Mathematical modeling analysis on the dynamics of university students animosity towards mathematics with optimal control theory

Animosity towards mathematics is a very common worldwide problem and it is usually caused by wrong information, low participation, low challenge tolerance, falling further behind, being unemployed, and avoiding the advanced math classes needed for success in many careers. In this study, we have considered and formulated the new SEATS compartmental mathematical model with optimal control theory to analyze the dynamics of university students’ animosity towards mathematics. We applied the next-generation matrix, Ruth-Hurwitz criteria, Lyapunov function, and Volterra-Lyapunov stable matrices to show local and global stability of equilibrium points of the model respectively. The study demonstrated that the animosity-free equilibrium point is both locally and globally asymptotically stable whenever the model basic reproduction number is less than unity, whereas the animosity-dominance equilibrium point is both locally and globally asymptotically stable when the model basic reproduction number is greater than unity. Finally, we applied numerical ode45 solvers using the Runge–Kutta method and we have carried out numerical simulations and shown that applying both prevention and treatment controls is the best strategy to minimize and possibly eradicate the animosity-infection in the community under consideration.

abilities are likely to have a more negative attitude towards mathematics and will not have the tendency to amend their mathematics skill. Researches indicated that poor attitudes towards mathematics are related to lower levels of achievement in the subject 26 . The main cause of the animosity towards mathematics especially in developing nations are attitudes of the learners, uninteresting lessons and low motivation of teachers attitude, confidence of learners, lack of teaching experiences, economic conditions, parents' educational level, teacher competency in math education, teacher's personality, intellectual factor, communication, stress, the pressure to perform well, lack of appropriate teaching methods, many numbers of student in the class over demanding tasks 21,27 . Even though, decreasing numbers of students choosing to study mathematics and natural science have been major issue, mathematics has been as a fundamental forerunner to success in worldwide modern society 5,21 . Since 1995, many theoretical researches have been conducted a series of international assessments of educational achievement in mathematics and science 31 .
Mathematical modeling approach using deterministic method 7,13,16,19,23,30 or stochastic method 11,18,[32][33][34] or fractional order method 20,24,25 have been scientific efforts to link and discover real world situations using mathematical models which have fundamental decision-making tools for the analysis of dynamics of communicable diseases and can be used for analyzing a number of real-world physical dynamical situations. The main purpose of this study is going to develop and analyze the new compartmental mathematical model with optimal control theory to analyze the animosity of higher institutions especially university students towards mathematics with tutorial and psychological treatment prevention and control measures. This mathematical model study is a first attempt in the thematic area, to carry out the proposed study we have faced a lack of mathematical modeling analysis literatures on higher institutions (university) students' animosity against mathematics and hence we have reviewed other literatures that are nearly related in our study. Mamo, Dejen Ketema, 2020 19 constructed a new Susceptible, Exposed, Infected Deny (SEID) racism expansion compartmental mathematical model, that explains the racism dissemination throughout a community under consideration. His study was mainly concentrated on the racism transmission on societal networks. He has been constructed and examined his proposed model theoretically and verified numerically. Applying the next-generation matrix approach, he has got the basic reproduction number of the model and the model basic reproduction number is almost correlated to the expansion of racism. Kooken, Janice W, et al. 20121 15 developed and presented the results of growth and establishment of the Cyclical Self-Regulated Learning (SRL) model simulation, a model of student knowledge and metacognitive learning mathematics experiences within an intelligent tutoring system (ITS). Their analysis results provide the establishment of the Cyclical SRL Model, emotion, performance in the ITS, and confirming the interplay of grit. Their simulation model enables mathematical simulations depicting a variety of student scope types and prevention styles and confirming deeper future student learning explorations. Van der Merwe, A, et al. 2018 35 presented interconnected algorithmic results that utilize mathematical programming models to bring forth and provide learning feedback in the form of academic performance status reports. Yadav, Anuradha, Prashant K. Srivastava, and Anuj Kumar 36 proposed and analyzed (PSQ) model where P stands potential smokers, S stands smokers and Q stands quitters for understanding of the dynamics of smoking behavior in a population. They first assumed that smokers are quitting smoking which is influenced by the level of determination of individuals. Higher degree of determination will lead to less chances of relapsing. Further the impact of education on potential smokers is also considered. They have shown that when individuals are educated about the fatality of the diseases caused by smoking they will refrain from smoking in future. Khyar, Omar, Jaouad Danane, and Karam Allali 14 explored mathematically the dynamics of giving up smoking behavior. For this purpose, they performed a mathematical analysis of a smoking model and suggested some conditions to control this serious burden on public health. The model under consideration describes the interaction between the potential smokers (P), the occasional smokers (L), the chain smokers (S), the temporarily quit smokers ( Q T ), and the permanently quit smokers ( Q P ). Existence, positivity, and boundedness of the proposed problem solutions are proved. Local stability of the equilibriums is established by using Routh-Hurwitz conditions. Moreover, the global stability of the same equilibriums is fulfilled through using suitable Lyapunov functional. In order to study the optimal control of their problem, they took into account a two controls' strategies, the government prohibition of smoking in public areas which reduces the contact between nonsmokers and smokers, the educational campaigns and the increase of cigarette cost which prevents occasional smokers from becoming chain smokers. The existence of the optimal control pair is discussed, and by using Pontryagin's minimum principle, these two optimal controls are characterized. Finally, numerical simulations are performed in order to check the equilibriums stability, confirm the theoretical findings, and show the role of optimal strategy in controlling the smoking severity. Alkhudhari, Zainab, Sarah Al-Sheikh, and Salma Al-Tuwairqi 2 derived and analyzed a mathematical model of smoking in which the population is divided into four classes: potential smokers, smokers, temporary quitters, and permanent quitters. In their model they studied the effect of smokers on temporary quitters. Two equilibriums of the model are found: one of them is the smoking-free equilibrium and the other corresponds to the presence of smoking. They examined the local and global stability of both equilibriums and support their results by using numerical simulations.
The main limitations of this study are: using non-modelling studies scholars shown us the problem is common in the community, however, in this study we faced problems to find more related literatures to the study and well organized real data of university students' animosity towards mathematics. Studies of researchers we reviewed did not considered and developed university students' animosity against mathematics mathematical modelling analysis and hence it makes this study has never been done by other scholars abroad. In our study we formulated and analyzed a new compartmental mathematical model with optimal control theory to minimize and possibly eradicate university students' animosity towards mathematics from the community and which shows the novelty of the study in the thematic area. Therefore, we are motivated by the limitations of the study to undertake this study and to fulfill the gap. The remaining part of the study is organized as; the deterministic model is constructed in Sect. Deterministic model formulation and is analyzed in Sect. Qualitative investigation of the deterministic model; the stochastic touch is explained in Sect. The transition probability; optimal control analysis is carried

Deterministic model formulation
In this study we considered the total number of students N(t) in a given time t and we divide it in to four disjoint classes. Those are susceptible, exposed, animosity infected and treated classes denoted by S, E, A, and T respectively. The state variables are describe as follows.
i. Animosity-susceptible students. Those are a group of student who are entering in to the university without a strong feeling of dislike or hatred of mathematics includes students who are interested in mathematics and became animosity infected whenever they meet with infected groups and it is denoted by S(t). ii. Animosity-exposed students. Those are a group of student who are taking mathematics course in the university and have closed contact with animosity student and hearing bad attitude towards mathematics and it is denoted by E(t). Those animosity exposed individuals may be animosity-infected or not. iii. Animosity-infected students. Those are a group of students who have a strong feeling of dislike or hatred or enmity that tends to display itself in mathematics and score lower grade in the university and it is denoted by A(t). iv. Animosity-treated students. Those are a group of students who are taking different tutorial and psychological treatments by their advisor and lecturers, we call it treated and denoted by T(t).

Basic assumptions and parameters definitions of the model
• The animosity-susceptible group S(t) increases by student who are newly entering in to the university with the rate , and by the number of student who are purely accept psychological treatment of their mentor and take their tutorial effectively with the rate of α , decreases due to the contact of animosity infected students with the rate β where the total number of students is constant with equal birth and death rate µ. • The animosity-exposed group E(t) increases by the students who have closed contact with animosity infected students and hearing about the animosity of mathematics with the rate β and decreases by the rate γ who are animosted mathematics. • The animosity infected group A(t) is increases by the rate γ and decreases by the tutorial and psychological treatment rate δ at any time t. • The treated group T(t) increases by the psychological treatment rate δ and decreases by the conversion rate α from treated class at any time t. • Students in all compartments are decrease by the natural death rate µ . Assume the total number of students is constant. • The parameter � = µN is the recruitment rate of freshman university students. Based on the model assumptions and descriptions above the flow chart (schematic diagram) of the flow of students is given in Fig. 1.
Using the model assumptions, parameter definitions and Fig. 1 the dynamical system of the study is given by www.nature.com/scientificreports/ We considered all model parameter values are non-negative and to analyze the mathematical model simply, the variables of the model (1) can be scaled as After substituting the time derivative of system (2) in system (1), we have got the new simplified system given by The total number of university students is constant, i.e.
Due to the scaled model of system (1), we have got the total number of students as Then the mathematically and biologically feasible domain of the system (3) is The model is both mathematically and biologically meaningful if following important lemma holds.

Lemma 1 The set is positively invariant to the system (3).
Proof Symbolize y(t) = (s(t), e(t), v(t), z(t)) T and the system (3) can be rewritten as.
Obviously the feasible region is closed set, and we want to prove if the initial condition f (0) ∈ � , then the solution f (t) ∈ � for all t ≥ 0.
If the dot product of f (y) the and normal vectors ( w 1 ;w 2 ;w 3 ;w 4 ; w 5 ) of the boundary lines are less than zero then y(t) ∈ � for all t ≥ 0.i.e.; (1)

Animosity-Dominance equilibrium point.
In the presence of animosity towards mathematics, the time dependent solution of the system (33) is said to be animosity-dominance equilibrium point denoted by E * given by E * = (s * , e * , v * , z * ) where after some steps of calculations we have got Suppose the characteristic polynomial of the matrix A is given by P A ( ) = det( I − A) = n + a 1 n−1 + . . . .. + a n and define k matrices as follows:

Then all eigenvalues have negative real parts, that is, the steady-state is stable if and only the determinants of all
Hurwitz matrices are positive etH j > 0 j = 1, 2, 3, ..k.

To show the local stability use linearization principle and Lyapunov function for the global stability of the system (3) equilibrium points respectively. Linearity of the system (3) determined by the help of Jacobean matrix is given by
The animosity-free equilibrium point is locally asymptotically stable if R 0 < 1.

Proof
The Jacobean matrix of the system (3) at the animosity-free equilibrium point is.
Proof The Jacobean matrix of the system (3) at the animosity-dominance equilibrium E * point is.
Then the eigenvalue of the system (5) is obtained from the characteristics equation The characteristics polynomial after simple simplification is where By Routh-Hurwitz stability criteria Theorem 3.1 above or by referring 13,23 , all eigenvalues of (6) have negative real parts see in 17,29 whenever R 0 > 1 . Hence whenever R 0 > 1, then the animosity-dominance equilibrium point is locally asymptotically stable. www.nature.com/scientificreports/ Thus dl dt < 0 , if R 0 < 1 and the equality dv dt = 0 holds if v = 0 and hence according LaSalle's invariant principle used in 30 the model animosity-free equilibrium point E 0 is globally asymptotically stable when R 0 < 1.

Theorem 3.4 The animosity-free equilibrium point is globally stable if
Global stability analysis of animosity-dominance equilibrium. The aim of this section is to investigate the global stability of E * in the positively invariant set with the aid of Volterra-Lyapunov stable matrices. To organize this, we express the Lyapunov function as: To establish the global stability of the animosity-dominance equilibrium point E * , we investigate the matrix A defined in Eq. (9) is Volterra-Lyapunov stable. We concisely explain the following basic definitions related to Volterra-Lyapunov stable matrices stated in 6 .
Let A nxn be a real matrices then. i. a nn > 0, Note: Ã is the (n − 1)x(n − 1) matrix, obtained by deleting the last row and column of the matrix A nxn .

Theorem 3.5
The square matrix A defined in (9) is Volterra-Lyapunov stable.
Proof Clearly A 44 > 0 , and the matrix when we delete last row and column of matrix −A is defined as: Based on Lemma 3, we state and prove that = −Ã , and C = − ∼ A −1 are diagonal stable in the following conditions to satisfied Lemma 3. Hence the matrix A is Volterra-Lyapunov stable.

Condition 1
The square matrix M defined in (10) is diagonal stable.

Proof
The diagonal stability of M , is guaranteed by the following steps: Step 1. It is clear that −Ã 33 = (µ + δ) > 0.
Step 2. By using Lemma 2, we need to show M is diagonal stable. From Eq. (10), we obtain Step 3. We show that It is easy to show that.
Hence ∼ M −1 is diagonal stable. Therefore in this condition we have proved that M = −Ã is diagonal stable.
Hence C is diagonal stable. Therefore from condition 2 we have proved C = − ∼ A −1 is diagonal stable. Finally from condition 1 and 2 we conclude M = −Ã and C = − ∼ A −1 are diagonal stable. Then finally due to the above basic definition, Lemmas and conditions, we have the following conclusions for the globally stability of the animosity-dominance equilibrium.

Stochastic touch for SEATS Model (1)
The transition probability. The transition probability is the probability of a stochastic process that transfer from state 1 to state 2. Here, the stochastic SEATS model consists of three random variables, i.e.
Prob (k1,k2,k3),))(s,e,a) (t, t + �t) = Prob{S(t + �t) = k 1 , www.nature.com/scientificreports/ Animosity outbreak probability. Animosity outbreak issue occurs when the number of animosity infected students' increases through time. The basic reproduction number (R 0 ) and the expected number of animosity infected students (n) are used as a criterion for the occurrence of animosity outbreaks in the long term. The basic reproduction number (R 0 ) is the number of susceptible students getting animosity infected after an animosity infected student is introduced into the group. The basic reproduction number (R 0 ) has the same definition as the expected number of infected students (n) , which is calculated using a probability method. In deterministic models, animosity outbreak issues occur when R 0 > 1 , whereas in stochastic models this is associated with a condition when the expected number of animosity infected students (n) > 1 . The two quantities (R 0 ) and n are obtained with different methods. The benchmark for the animosity outbreak probability is determined by the expected number of animosity infected individuals (n) , not from the basic reproduction number (R 0 ) . Based on this process, the SEAS model has the following disease extinction probabilities: and as a result, the probability of animosity outbreak occurring is.
Finally, in this study we considered large number of students for numerical simulation; for the remaining analysis part we have used the deterministic model rather than the stochastic touch.

Optimal control analysis of the deterministic model
In this section, we formulated and analyzed a mathematical model with optimal control method 34 to identify the best control strategy that reduces the number of animosity-infected individuals with minim total cost used during the interventions. The objective is to find the optimal values u * = u * 1 , u * 2 of the controls u = (u 1 , u 2 ) such that the associated state trajectories E * = (S * , E * , A * , T * ) are solution of the system (1) in the intervention time interval 0, T f with given initial conditions and minimize the objective functional. The control u 1 (t) represents the efforts on preventing animosity infection that helps to reduce contact rate of animosity, u 2 (t) is the control relates towards improvement of animosity recovery period. This is equated towards implementation of proper treatment and psychological counseling policies for animosity towards mathematics, so that recovery period can be improved by providing right treatment and counseling such that 0 ≤ u 2 (t) ≤ 1.

Then the system (1) is changed into
With the corresponding initial conditions For this, our optimal control problem is to minimize the objective functional where I(S, E, A, T, u) = w 1 A + B 1 2 u 2 1 + B 2 2 u 2 2 , measures the current cost at time t. The coefficient w 1 is positive weight constant that characterizes the cost associated with minimizing the animosity infected individuals and B 1 2 and B 2 2 are the measures of relative costs of interventions associated with the controls u 1 and u 2 , respectively, and also balances the units of integrand. In the cost functional, the term w 1 A refers the cost related to animosity infected class.
The set of admissible control functions is defined by More precisely, we seek an optimal control pair Characterization of the optimal control. In this section, we present optimality conditions for the optimal control problem defined above and detail its properties. According to Pontryagin's Maximum Principle in 4,13,37 , if u * (t) ∈ � u is optimal for dynamical system (10) with initial value (11) and (14) with fixed final time T f , then there exists a non-trivial absolutely continuous mapping : 0, T f → R 4 , = ( 1 (t), 2 (t), 3 (t), 4 (t)) called the adjoint vector, such that.
(1) The Hamiltonian function is defined as also holds true. In the next result, we have discussed characterization of optimal controls and adjoint variables.
Theorem 5.1 Let u * = u * 1 , u * 2 be the optimal control and (S * (·), E * (·), A * (·), T * (·)) be the associated unique optimal solutions of the optimal control problem (10) with initial condition (11) and objective functional (12) with fixed final time T f (13). Then there exists adjoint function * i (·), i = 1, ..., 4 satisfying the following canonical equations: with transiversality conditions Moreover, the corresponding optimal controls u * 1 (t) and u * 2 (t) are given by From the previous analysis, to get the optimal point, we have to solve the system

Numerical analysis
Numerical methods. In this study we have applied the numerical MATLAB ode45 programming code written by using the Runge-Kutta methods which are generally more powerful method for solving non-stiff ordinary differential equations.
Numerical simulations for the deterministic model. In this section, we convey the numerical simulations to verify the theoretical results of our mathematical model (3). Particularly, some numerical justifications are considered to illustrate the theoretical analysis and results of the preceding sections. Here we assume the parameter values for numerical simulations that are not from real data, since there is the lack of mathematical modeling analysis literatures which have been done to study the dynamics of animosity towards mathematics.
The initial values of the model (2) are positive, i.e.,s(0) > 0, e(0) ≥ 0, v(0) ≥ 0, z(0) ≥ 0 since the model represents human population. To understand the dynamics of students' animosity towards mathematics, we need to assume parameter values and extremely analyze the model (3), and further explore how these parameters influence on the spread of animosity. For instance parameter describes the possibility that s(t) converts to be an e(t) due to the presence of single i(t) per unit time. It is essentially determined by the rate of an i(t) spreading animosity.
Numerical simulation when the threshold quantity R 0 less than unity. In this subsection let us consider the mathematical model (3)  Here from Fig. 2 we can justify that the model animosity-free equilibrium point is both locally and globally stable whenever R 0 = 0.97 < 1 , practically it means that animosity towards mathematics eradicate from the group of students in the near future (after 150 months) through the community.
Numerical simulation when the threshold quantity R 0 greater than unity. In this sub-section let us consider the parameter values as = 0.5 , γ = 0.04 , α = 0.03 , µ = 0.02 , δ = 0.01 , then we have got the basic reproduction number of model (2) as R 0 = 11.1 , and we confirmed that the animosity-dominance equilibrium point is locally and globally stable.
Here from Fig. 3 we can justify that the model animosity-dominance equilibrium point is both locally and globally stable whenever R 0 = 11.1 > 1 , practically it means that animosity of students towards mathematics exists uniformly after 150 months in the group of students throughout the community.
Effect of the transmission β on the animosity-exposed students. Here we perform numerical simulation of animosity-exposed students "e" with parameters values given by γ = 0.035 , α = 0.03, µ = 0.02, δ = 0.01 , and variable transmission rate β . Here from Fig. 4 we can justify that the animosity-exposed number of students "e" is going up whenever the contact rare or transmission rate β increases.   www.nature.com/scientificreports/ offs the animosity infected students and the positive constants B 1 , and B 2 represent the weight constants for the efforts on preventing animosity infection and improved animosity treatment respectively. The values assigned to the weight constants are w 1 = 1 , B 1 = 100 , and B 1 = 100 as given in 3 .
Optimal control simulation when both controls are applied. Here we carried out numerical simulation on the model (1) animosity-infectious class ( A ) with parameters values given by β = 0.04 , γ = 0.035 , µ = 0.02 , δ = 0.01 by applying both controls (i.e. u 1 = 0 and u 2 = 0 ), it results in significant difference in the number of infections under the application of controls which is shown in Fig. 7. We note that for the class of animosity infected students there is huge decrease in the number of infections when both controls are applied, and the number nears to zero after 250 months.
Optimal control simulation when exclusively one of the controls is applied. Here we carried out numerical simulation on the model (1) animosity-infectious class ( A ) with parameters values given by β = 0.04 , γ = 0.035 , µ = 0.02 , and δ = 0.01 , by applying only treatment control (i.e. u 1 = 0 and u 2 = 0 ), it results in significant difference in the number of infections under the application of treatment only control which is shown in Fig. 8. We note that for the class of animosity infected students there is a decrease in the number of infections when treatment only control is applied, and the number nears to zero after 150 months and whenever without also  www.nature.com/scientificreports/ treatment control i.e., when we do not apply both controls measures the animosity infectious class is increasing and becomes constant after 150 months.
Here we carried out numerical simulation on the model (1) animosity-infectious class ( A ) with parameters values given by β = 0.04 , γ = 0.035 , µ = 0.02 , δ = 0.01 , by applying only prevention control (i.e. u 1 = 0 and u 2 = 0 ), it results in significant difference in the number of infections under the application of prevention only control which is shown in Fig. 9. We note that for the class of animosity infected students there is a decrease in the number of infections when prevention only control is applied, and the number nears to zero after 200 months and whenever without also prevention control the animosity infectious class is increasing and becomes constant after 200 months. Finally, parameter values used to simulate Figs. 7, 8, and 9 have a fundamental impact to show how to control the students' animosity towards mathematics using suitable control measures applied in the model formulation.

Conclusion
In this study we performed a detailed study on the deterministic with some stochastic touch mathematical model on the students' animosity towards mathematics with optimal control approach comprising of four mutually exclusive compartments. The study began with a detailed analysis covering equilibriums, basic reproduction number, and stability analysis of the disease-free and endemic equilibriums. From the analysis on students animosity towards mathematics model we have proved that the disease-free equilibrium and the endemic www.nature.com/scientificreports/ equilibrium are locally asymptotically stable when the basic reproduction number ( R 0 ) is less than unity and greater than unity respectively. The theoretical results are then verified by optimal control analysis and numerical simulations respectively and results shown us the following: • The model solutions converging to the animosity infection-free equilibrium point whenever R 0 = 0.97 < 1, biologically the animosity infection will be eradicate in the near future. • The model solutions converging to the animosity-dominance equilibrium point around 250 months whenever R 0 = 11.1 > 1, biologically the animosity infection will be spreading in students. • The animosity-exposed number of students "e" is going up whenever the contact rare or transmission rate β increases, whenever the progression rate γ increases the number of animosity-infected students " v " going up, and whenever the value of treatment rate δ increases then the number of animosity-infected students " v " going down.
Optimal control analysis performed by including two control parameters, one is the prevention mechanism (u 1 (t)) and the other is treatment of animosity infection mechanisms such as psychological treatment, tutorial methods (u 2 (t)). A significant difference in the number of infections under the application of controls was shown in the results discussions. The effect of these controls individually was also performed and when we applied both controls together, the reduction in the number of infected population was quite large. Hence, this analysis focused on the need of improved prevention mechanisms and treatment to minimize and possibly to eradicate the animosity from students population. Hence, the results from the study imply acceleration in the animosity spread if proper prevention and treatment are not implemented. On a concluding note, the study further suggests that in the times of animosity spreading prevention and treatment should not be neglected. Although some literatures have acknowledged university students animosity towards mathematics, no study has attempted to explore the idea with mathematical modelling approach. This study attempts to contribute to the educational research by formulating and analyzing a mathematical modelling approach with optimal control approach on animosity towards mathematics which allows for deep understanding and a more comprehensive idea on the thematic area and it is distinguished from previous studies by using mathematical modelling approach as infectious disease to examine animosity towards mathematics.